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^ ; ABSTRACT 

■ We present the first AGN census in a sample of 6f galaxies selected at 70 /im, a wave- 
On' length which should strongly favour the detection of star-forming systems. For the pur- 
q , pose of this study we take advantage of deep Chandra X-ray and Spitzer infrared (3.6- 
5^ ■ 160 /im) data, as well as optical spectroscopy and photometry from the Deep Extra- 

' galactic Evolutionary Probe 2 (DEEP 2) survey for the Extended Groth Strip (EGS) 

Ci . field. We investigate spectral line diagnostics ([OIIIJ/H/3 and [NeIII]/[OII] ratios, H5 

Balmer absorption line equivalent widths and the strength of the 4000 A break), X- 
| ray luminosities and spectral energy distributions (SEDs). We find that the 70 /im 

■ sources are undergoing starburst episodes and are therefore characterised by a pre- 
00 ' dominance of young stars. In addition, 13 per cent of the sources show AGN signatures 
£f~) , and hence potentially host an AGN. When the sample is split into starbursts (SBs, 
O ■ 10 10 <L IR <10 n L Q ), Luminous InfraRed Galaxies (LIRGs, 10 11 < L m < 10 12 L ) 
m ; and UltraLuminous InfraRed Galaxies (ULIRGs, 10 12 < h IR < 10 13 L ), the AGN 

. fraction becomes 0, 11 and 23 per cent respectively, showing an increase with to- 

— — ■ tal infrared luminosity. However, by examining the sources' panchromatic SEDs, we 

| conclude that although the AGN is energetically important in 1 out of 61 objects, 

. all 70 /xm-selected galaxies are primarily powered by star-formation. When compared 

to a sample of DEEP 2 galaxies in the same redshift range and with similar optical 
colours, we find that the 70 /im population is characterised by younger stellar ages and 
a higher AGN incidence, indicating that strongly star-forming populations might be 
key in studying the relationship between black hole and stellar growth. 

a ' 

1 INTRODUCTION 1997; Appleton et al. 2002; King 2005; Springel, Di Mat- 

teo & Hernquist 2005) . Although, disentangling the relative 

Understanding the relationship and relative contribution of starburst/AGN contributions to the infrared energy budget 

the fundamental energy sources in the Universe, AGN ac- is key for identifying the processes that drive galaxy evo- 

cretion and star formation, is far from trivial. Data from lution, such a task has proved extremely challenging and 

three generations of infrared observatories, IRAS, ISO and any substantial progress has been slowed by the difficulty 

Spitzer, have revealed that the most powerful activity in the of identifying obscured AGN at high redshift. The presence 

Universe is dust-enshrouded. Obscured galaxies are respon- of dust in the host galaxy complicates matters further, as 

sible for roughly half the cosmic energy density originat- the relative contribution of stars and AGN to the bolometric 

ing from stars and AGN: the Cosmic InfraRed Background luminosity becomes an elusive quantity. 

(CIRB). Observational and theoretical results have placed The main contributors to the CIRB, Luminous and Ul- 

the history of star formation and black hole accretion in traLuminous Infrared Galaxies (LIRGs and ULIRGs), are 

parallel evolutionary paths, linked to the close relationship characterised by large amounts of dust, emission from which 

between central black holes and galaxy properties (Magor- completely dominates their energetic output. As dust grain 

rian et al. 1998), a key aspect of identifying the processes formation and destruction is a dynamical, stellar-related 

that drive galaxy evolution (e.g. Sanders et al. 1988; Nor- process with the lifetime of a typical grain of the order 

man & Scoville 1988; Rowan- Robinson & Crawford 1989; of a few hundred Myr (Draine 2003), the presence of large 

Rowan- Robinson 1995), with speculations of a complex two- amounts of dust combined with a high infrared output im- 

way feedback method whereby AGN and starburst activ- plies rapid, ongoing star- formation and/or intense starburst 

ity is inter-regulated (e.g. Umemura, Fukue & Mineshige episodes. With data from the Multiband Imaging Photome- 
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ter for Spitzer (MIPS), the properties of these populations 
have been examined out to z ~ 2, however it is only with 
the far-IR (70 and 160 fim) bands that greater relevance to 
star-formation can be achieved, as they probe nearer the 
SED peak of a star-forming galaxy (~ 40-150 /im). Because 
of its higher sensitivity, the MIPS 70 /jm band has been more 
widely exploited than the 160 pm band, with 70 /jm popu- 
lations found to predominantly consist of infrared-luminous 
(L/i? > 10 10 ) and dust-rich objects, mainly in the 0.1<z<1.5 
redshift range (Symeonidis et al. 2008; 2009). As a result, 
setting the initial sample selection at 70 /im is ideal for 
studying the relationship (if any) between obscured star- 
formation and supermassive black hole growth, as well as 
evaluating the degree of coeval activity between AGN and 
star-formation during a galaxy's lifetime. 

Numerous studies have revealed that a large part of the 
bolometric luminosity in IR-luminous sources could origi- 
nate from AGN in some cases prevailing over the starburst 
output (e.g. Gregorich et al. 1995; Genzel et al. 1998; Tac- 
coni et al. 2002), with the AGN contribution seen to increase 
as a function of total infrared luminosity (e.g. Lutz et al. 
1998; Fadda et al. 2002; Brand et al. 2006). Identifying the 
presence of AGN in obscured galaxies requires the combi- 
nation of various disgnostics such as mid-IR colours (e.g. 
Lacy ct al. 2004; Stern et al. 2005; Barmby et al. 2006), 
SED shape (e.g. Alonso-Herrero et al. 2006), X-ray lumi- 
nosity (e.g. Pompilio, La Franca & Matt 2000) and radio 
luminosity (e.g. Donley et al. 2005). Optical spectrosopy is 
also a powerful tool (e.g. Osterbrock 1989; Tresse et al. 1996; 
Takata et al. 2006), with recent evidence pointing towards 
evolutionary differences, rather than solely orientation ef- 
fects (e.g. Antonucci 1993), being responsible for the vari- 
ety of AGN spectroscopic signatures (e.g. Hasinger 2008), 
resulting in a mixture of both type I and type II spectra. 
In addition to the traditional 'torus-obscured' AGN, stud- 
ies have revealed the potential of 'host-obscured' AGN as a 
plausible scenario, where light from the nucleus is also ob- 
scured by the dusty star-forming host galaxy (e.g. Martinez- 
Sansigre et al. 2006). For such sources, optical spectroscopy 
is a powerful diagnostic with a plethora of emission lines 
both from the AGN and host, providing important clues on 
the environments in which they originate (e.g. Brand et al. 
2007). 70 /im populations are at a considerable advantage 
for multiwavelength spectroscopic and photometric follow- 
up: the sensitivity of the MIPS 70 /im band is such that only 
the brightest sources are selected at each redshift and hence 
are very likely to have counterparts over a wide wavelength 
range. Moreover, the redshift distribution of 70 fxra popula- 
tions allows many key emission and absorption lines to be 
within a typical optical grating range. 

Due to the effects of dust obscuration, AGN censuses 
in infrared and sub-mm populations have mainly been ap- 
proached through mid-IR and X-ray surveys (e.g. Vig- 
nati et al. 1999; Almaini, Lawrence & Boyle 1999; Fabian, 
Wilman & Crawford 2002; Alexander et al. 2005; Martinez- 
Sansigre et al. 2005; Matute et al. 2006). X-ray surveys 
have resolved most of the 0.5-10 keV Cosmic X-Ray Back- 
ground (CXRB) and attributed it to emission from accreting 
super-massive black holes (e.g. Shanks et al. 1991), how- 
ever lack of sensitivity has not allowed the detection of 
the majority of extragalactic X-ray sources contributing at 
low (< 10 -16 erg/s/cm 2 ) fluxes. These are most likely high- 



redshift, star-forming and starburst galaxies, possibly con- 
taining low luminosity or deeply obscured AGN (e.g. Ptak et 
al. 2003; Franceschini et al. 2003; Bauer et al. 2004; Brandt 
& Hasinger 2005). A large obscured AGN population has 
long been predicted both from statistical analysis and num- 
ber counts of AGN samples as well as from CXRB models 
(e.g. Comastri et al. 1995; Worsley et al. 2005), especially 
since the steep broad-line QSO spectrum cannot explain the 
intensity of the hard X-ray background (e.g. Fabian et al. 

1998) . This implies that most accretion-generated energy 
density in the Universe takes place in obscured objects, mak- 
ing obscured AGN the main candidates for the origin of the 
hard CXRB, which peaks at ~30keV and, as yet, has not 
been resolved (e.g. Setti & Woltjer 1989; Fabian & Iwasawa 

1999) . 

The aim of this paper is to examine the AGN content 
of a sample of IR-luminous 70 /im-selected galaxies from the 
Extended Groth Strip (EGS) field. We combine Chandra X- 
ray data and optical spectroscopy in order to constrain the 
AGN fraction and identify the dominant process responsible 
for the sources' total energy budget. Section 2 introduces the 
sample and gives an analysis of IR/optical colours. In section 
3 we examine the sample's X-ray properties, whereas in sec- 
tion 4 the analysis focuses on optical spectra and emission 
line ratios. The AGN contribution is discussed in section 5 
and our summary and conclusions are presented in section 
6. Throughout we employ H = 70 kms^Mpc -1 , Q M = 0.3 
and tt A = 0.7 (Spergel et al. 2003). 



2 THE SAMPLE 

2.1 Initial selection and previous work on the 
70 /im sample 

This work is based on Guaranteed Time Observations 
(GTO) of the Extended Groth Strip (EGS) field (~0.5deg 2 ) 
(Davis et al. 2007) by Spitzer's (Werner et al. 2004) far- 
IR photometer MIPS (Rieke et al. 2004). The initial selec- 
tion was made at 70 /jm, where 178 sources were retrieved 
down to ^4mJy (5a), with photometric completeness at 
~ lOmJy. We focus on a subset of this population, the 114 
sources in the overlap areas of the MIPS 24, 70, 160 /jm 
and InfraRed Array Camera (IRAC) 8 /jm images and fur- 
ther narrow our working sample to 61 sources with optical 
photometry and spectroscopic redshifts from the Deep Ex- 
tragalactic Evolutionary Probe 2 (DEEP 2) survey (see sec- 
tion 2.2). The MIPS 70 pm catalog was cross-correlated to 
the MIPS 24 /im and IRAC catalogs in two steps: first us- 
ing a ~ 10 "radius for the MIPS70/MIPS24 matching and 
subsequently using a ~ 2.5 "radius for the MIPS24/IRAC 
matching. For data reduction and source extraction we refer 
the reader to Symeonidis et al. (2007; 2008; 2009, hereafter 
S07; S08; S09), where the infrared properties of the 70 /jm 
sample were also examined in detail. A succinct summary 
of our results from previous work follows: 

In S08 and S09, we fitted the available IR photometry 
(8-160 /im) with the Siebenmorgen & Krugel (2007) tem- 
plates and obtained estimates for the total infrared lumi- 
nosity (Lir) in the 8-1000 //m range. We found character- 
istic luminosities of IR-luminous galaxies between 10 10 and 
10 14 L Q , of which ~11% are Starbursts (SBs, 10 10 < L IR 
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Figure 1. The R-band magnitude distribution for the 109 70 /im 
sources — solid line: the 47 sources without redshifts; dashed 
line: the 61 objects with DEEP 2 redshifts. The exclusion of the 
brightest sources and spectroscopic weighting scheme, cause dis- 
crepancies in the two distributions at the faint and bright ends. 



< 10 1 ), ~62% are Luminous InfraRed Galaxies (LIRGs, 
10 11 < L IR < 10 12 ) and ~26% are UltraLuminous InfraRed 
Galaxies (ULIRGs, 10 12 < L IR < 10 13 ). We also identi- 
fied one HyperLuminous InfraRed Galaxy (HyLIRG, Lir 
> 10 13 ) at z=1.9, which we remove from subsequent analy- 
sis as it is unique and will be treated separately in a later 
paper. The optical/IR SEDs for all, but one, objects are 
starburst-type, with a strong optical/near-IR stellar bump, 
an inflection in the near-IR and elevated flux in the infrared. 
One object, identified with red near-IR colours in all IRAC 
colour bands, has an SED that shows an elevated near-IR 
power-law continuum instead of an inflection in that region 
(see sections 5 and 5.1). Infrared and radio derived star for- 
mation rates (SFRs) were calculated in S07 and were found 
to be in the ~ 5-1000 Mq/yi range, with the majority of 
the objects having SFR > 50 M0/yr. SFRs calculated using 
optical emission lines not corrected for dust extinction, on 
average fell short by a factor of 50, implying that the opti- 
cal part of the spectrum is not an accurate tracer of star- 
formation in these systems, especially since in addition to 
high extinction, IR-luminous sources will likely host deeply 
embedded star-forming regions that do not contribute at all 
to the optical energy budget. 



2.2 Optical selection 

The EGS benefits from extensive observational coverage 
via the All- Wavelength Extended Groth Strip International 
Survey (AEGIS), which provides optical spectroscopy and 
multi- wavelength photometry over an area of ~ 1 deg 2 (see 
Davis et al. 2007 for a summary of the AEGIS photomet- 
ric datasets). From the 114 objects with infrared photome- 
try, we select 109 with area coverage in the B, R, I bands, 
from the DEEP 2 Survey (Coil et al. 2006; Faber et al. 
in prep.). Taking advantage of IRAC's high positional ac- 
curacy (~0.l"), we cross-correlate the 70 (jm sample with 
DEEP 2 using the IRAC coordinates for our sources and 
an appropriate matching radius of < 1.5 "corresponding to 
IRAC's pixel size. Currently, ~ 75 per cent of the MIPS 



Figure 2. R-band magnitude versus redshift. 
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Figure 3. 70 fim flux density versus redshift. 



70 /im strip is covered by DEEP 2 and this will increase to 
100 per cent with the DEEP 3 cycle of observations. Where 
possible, we also supplement our dataset with photometry 
from the Galaxy Evolution Explorer (GALEX), the Palo- 
mar Wide-field InfraRed Camera (WIRC) and the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) in the 
FUV, NUV, g, r, i, z, J and K bands and Spitzer's IRAC in 
the 3.6, 4.5 and 5.8 /tm bands. 

Apart from the exclusion of stars, the photometric part 
of DEEP 2 did not pose any constraints on the EGS. On 
the other hand, the spectroscopic part of DEEP 2 targeted 
sources only within the 18.5 <R< 24.1 (AB) magnitude 
range. Although the retrieval of redshifts was random at ~60 
per cent sampling, a weighting scheme was applied, giving 
lower weights to z<0.75 objects in order to sample a wide 
range of luminosities and roughly equal numbers of galaxies 
below and above z=0.75; see Faber et al. (2007); Davis et 
al. (2007). For a reliable redshift, two significant features, 
such as the [Oil] doublet, were required in order to provide 
a match and hence most redshift failures were from high- 
redshift objects (z > 1.4) which did not have strong features 
in the DEEP 2 spectral range. At the time of writing there 
are 62 good quality redshifts and spectra available for the 
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Figure 4. U-B colours vs absolute B magnitude, Mb (AB), with the 70 fim sample in black open triangles, blue filled circles and red 
asterisks(SBs, LIRGs and ULIRGs, respectively) and the control sample as green filled squares. All ~ 11,000 EGS source from the DEEP 2 
spectroscopic survey are shown in the background as black points. 



70 fim EGS population, in the 0.17<z< 1.9 range, with a 
mean of 0.57 and a median of 0.48. Note that we remove the 
HyLIRG from the redshift sample (see section 2.1) with the 
final number of sources being 61. Hereafter the term '70 fim 
sample' refers only to the sources with redshifts. 

In figure 1 we compare the R-band magnitude dis- 
tribution of the redshift sample to the sources with no 
redshifts. Although there is high overall consistency, a 
few differences are evident originating from the spectro- 
scopic selection criteria outlined earlier. Firstly, there is an 
under-representation of bright sources: there are 10 R<18.5 
sources in total and, at 60 per cent redshift sampling, ~ 6 
of those should have made it in the redshift sample. The 
R-band - redshift relation (Fig. 2) shows that these are at 
low redshift. Secondly, on the faint end, there are 4 objects 
fainter than R=24.1 in the full sample, i.e. 2-3 should have 
also appeared in the redshift sample. The effects of the spec- 
troscopic selection are further examined in figure 3, where 
the 70 /J,m flux density (fro) is plotted as a function of red- 
shift. In contrast to the clear R-z relation (Fig. 2), £70 does 
not show signs of a trend with redshift, implying that the 
DEEP 2 spectroscopic criteria do not have any affect on the 
infrared properties of the 70 /im sample. 

2.3 The Control Sample 

Due to the unique nature of the 70 ^m population (see sec- 
tion 2.1) and as this work focuses on determining the AGN 
fraction, it is more insightful to do so in relation to galaxies 




0.4 0.6 0.8 1.0 
spectroscopic redshift 

Figure 5. The redshift distribution of the 70 /im sample (red 
solid line) and the control sample (black dashed line) 



of comparable properties in a similar redshift range. With 
that in mind, we assemble a control sample of DEEP 2 EGS 
sources: for each 70 fj,m galaxy with a DEEP 2 redshift, we 
find up to three others with comparable redshifts, rest-frame 
B-band luminosities and U-B colours (see Willmer et al. 2006 
for a discussion of ^-corrections) — 'comparable' is quan- 
tified as Az ^0.03, AM S ^0.4 and A(U-B) sJO.l (see figures 
4 and 5). 

Figure 4 shows the variation of U-B colour versus Mb 
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(AB) — the Colour Magnitude Diagram (CMD) — for the 
control and 70 (im sample; the colour distribution for all 
EGS sources in the DEEP 2 spectroscopic survey is also 
shown in the background. The CMD represents galaxy evo- 
lution from the blue cloud — young stellar content — to 
the red sequence — evolved stellar populations — associ- 
ated with morphological changes from spirals to ellipticals 
and mediated by a transitional green valley stage (e.g. Bell 
et al. 2004) . The 70 /xm sample is bright in the optical but 
has red U-B colours; there is substantial contribution to the 
red sequence, but more sources seem to be in the upper blue 
cloud/green valley region. 

An alternative way to select the control sample could 
be by stellar mass. However, due to lack of complete 
optical/near-IR datasets for the EGS, mass information for 
the DEEP 2 galaxies (see Conselice et al. 2007) is only avail- 
able for about 64 per cent of the 70 /im population. The 
subset of 70 fim objects with available stellar masses are 
characterised by a narrower distribution of higher average 
mass (<log M> =10.87), compared to the control galaxies 
(<log M> =10.23) (C. Conselice, private communication). 
The stellar mass distribution of the 70 fim galaxies is consis- 
tent with previous studies of infrared galaxies (e.g. Caputi 
et al. 2006a; 2006b) and comparable to those of early-types 
in the z<l redshift range (e.g. Ferreras et al. 2009). For 
our purpose, colour/magnitude selection of the control sam- 
ple is more appropriate than stellar mass, because i) stellar 
masses are not available for the whole 70 /jm sample and ii) 
colour/magnitude selection avoids a control sample domi- 
nated by high mass evolved galaxies with little or no star- 
formation. 

Note that the optical luminosities and colours are not 
extinction-corrected, which might imply that the intrinsic 
properties of the 70 /j,m sample could be better matched to 
objects in the more optically luminous part of the CMD 
or those whose red colours are principally a consequence of 
dust-reddening rather than stellar age. Notwithstanding, the 
lack of correlation in Fig. 4 between U-B colour and infrared 
luminosity for the 70 /im sample, implies that subsequent 
comparisons of each IR luminosity class with the respective 
control galaxies will be fair, since the latter span a large part 
of the CMD and hence a large range in properties. 

2.4 X-ray Data 

The EGS was surveyed with the Chandra Advanced CCD 
Imaging Spectrometer (ACIS) to a depth of 200 ks, over an 
area of 0.67 deg 2 (Georgakakis et al. 2006; Davis et al. 2007; 
Laird et al. 2009). The AEGIS X-ray survey and final data 
products are made available in Laird et al. (2009, L09). L09 
describes the cross-matching between the DEEP 2 and X-ray 
AEGIS catalogs; we identify the X-ray counterparts of our 
70 /jm sample using DEEP 2 coordinates which were in turn 
matched to the IRAC coordinates of our sources (see section 
2.2). We use X-ray data in the soft (0.5-2 keV) and hard (2- 
lOkeV) bands (flux limit of 1.1 x 10~ 16 and 8.2 x 10~ 16 
erg/s/cm 2 respectively), choosing fluxes calculated using a 
Bayesian method, more reliable for faint sources as it cor- 
rects for the Eddington bias — for details on source extrac- 
tion and flux estimation see L09. 

Although the Chandra map covers the whole of the 
MIPS 70 /im strip, only about 16 per cent of the sources 



have individual X-ray detections (see table 4). Estimates 
for the hardness ratio, calculated as HR = |f^§, where H 
and S refer to the hard (2-7keV) and soft (0.5-2keV) count 
rates (counts/s/cm 2 ), are also shown. If there is only a de- 
tection in one band, HR is calculated using the 3u upper 
limits. From the 61 sources in the sample, only 10 have a 
confirmed detection in the full band (0.5-8 kev), 2 of which 
are not detected in any further bands, 4 have detections only 
in the soft band, 3 in both the soft and hard bands and 1 
in only the hard band. For the remaining objects we per- 
formed stacking (see table 1) , after splitting the sample into 
3 groups according to the total infrared luminosity measure- 
ments. The Tf-corrected X-ray luminosities are calculated by 
assuming a photon index of F= 1.9 and the stacked fluxes 
are converted to luminosity using the mean redshift of the 
group, 0.22 for SBs, 0.5 for LIRGs and 0.87 for ULIRGs. 



3 THE X-RAY/IR RELATION 

Figure 6 compares X-ray and infrared luminosities for the 
70 lira sample. Given that there are several sources of X-ray 
luminosity associated with stellar-related processes (e.g. X- 
ray binaries, supernova remnants, hot gas, starburst-driven 
outflows), X-ray and infrared emission are expected to be 
strongly correlated in star-forming galaxies (e.g. Fabbiano 
1989; Griffiths & Padovani 1990; Ranalli, Comastri & Setti 
2003; Franceschini et al. 2003), for which Lo.s-iofceV ranges 
between 10 40 -10 42 erg/s (Kim, Fabbiano & Trinchieri 1992a; 
1992b; Nandra et al. 2002; Laird et al. 2005). The bound- 
aries between star-formation and accretion typically extend 
over a couple of orders of magnitude, nonetheless a hard X- 
ray luminosity of L 2 _iofc e v > 10 42 erg/s is thought to be an 
indication that an AGN is the dominant emitter of X-rays 
(e.g. Georgakakis et al. 2007). 

We evaluate the properties of our sources against the 
empirical soft and hard X-ray/IR relations derived for IRAS 
galaxies in the 9 < log Leo < 11.5 range from Griffiths & 
Padovani (1990) and for local star- forming SBs and LIRGs 
from Ranalli, Comastri & Setti (2003, hereafter RCS03), as 
well as the hard X-ray/IR relation derived for a sample of 
local ULIRGs from Franceschini et al. (2003, hereafter F03). 
From figure 6, it is clear that the survey limits are displaced 
from the star-forming relations by up to 2 orders of magni- 
tude, suggesting that a 200 ks X-ray survey is not on average 
sensitive to distant star-forming galaxies and hence more 
likely to detect objects with a strong AGN contribution to 
their X-ray luminosity. 

In the soft X-rays (Fig. 6, top panel), the AGN/star- 
formation separation is not clear: there is large scatter in 
Lo.5-2koV, partly due to the fact that soft X-rays are subject 
to stronger attenuation by gas and dust, the degree of which 
is likely to vary from one system to another and as a result 
the contribution of the AGN to the soft X-rays may not be 
evident. In fact, even some of the ULIRGs in the F03 sam- 
ple whose X-ray emission is classified as AGN-dominated by 
F03 seem consistent with the star-formation relations. On 
the other hand, the 4 sources detected in the hard X-rays 
(objects 59, 67, 83 and 122, Fig. 6, lower panel) have lu- 
minosities of L2-iofeeV >10 42 erg/s, 2-3 orders of magnitude 
higher than expected from the IR-X-ray relation for star- 
forming systems, placing them into the AGN regime. Their 
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Table 1. Table of stacked fluxes (erg/s/cm 2 ) and hardness ratios for the 70 jim sample split into 
3 luminosity classes (starbursts, LIRGs and ULIRGs; see section 2.1 for details on the sources' 
infrared properties). The second column indicates how many objects were used in the stacking — 
some were removed either because they were too close to a source which could have contaminated 
the signal or they were associated with an X-ray source. The SB group does not have a formal 
detection in the hard band, so we quote the 3ct upper limit. The hardness ratio is derived using 
the counts in the soft band and 3cr upper limit counts in the hard band. 

Class Objects used fo.5-2 12— 10 HR 

Starbursts 8 out of 9 8.67E-17 (8cr) 2.51E-16 (3<r upper limit) -0.3 (upper limit) 
LIRGs 29 out of 35 5.04E-17 (9cr) 1.6E-16 (4cr) -0.2 

ULIRGs 9 out of 17 1.11E-16 (12ct) 1.84E-16 (3a) -0.5 
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Figure 6. Soft (top panel) and hard (lower panel) X-ray luminosity (erg/s) vs Ljr (erg/s lower x-axis and Lq upper x-axis) for the 70 fim 
sample. The colour scheme is green for SBs, blue for LIRGs and red for ULIRGs. Objects with a detection in at least one band have more 
reliable 3cr upper limits in other bands (large coloured symbols), whereas for the remaining sources the upper limits are calculated from 
the survey flux limits (see section 2.4, black symbols). Filled circles correspond to luminosities calculated using the stacked 0.5-2 keV and 
2-10 keV band flux and open squares are objects with X-ray detections. In the top panel, the lines are the star-forming relations from 
Griffiths and Padovani (1990) for star burst /interacting galaxies (dashed line) and IRAS galaxies (dot-dashed line) and from Ranalli, 
Comastri & Setti (2003) for HII galaxies (solid line). In the lower panel the lines indicate the star- forming relations from Franceschini 
et al. 2003 (dashed line) for a sample of ULIRGs and Ranalli, Comastri & Setti (2003) for HII galaxies (solid line). Overplotted is 
the ULIRG sample from Franceschini et al. (2003); ULIRGs which are SB-dominated in the X-rays (diamonds), ULIRGs which are 
AGN-dominated in the X-rays (triangles) and ULIRGs which are SB/AGN hybrids in the X-rays (asterisks). The open squares are the 
HII galaxies from Ranalli, Comastri & Setti (2003). The small cross is Arp220 and the large vertical cross is NGC6240, where the extent 
of the cross represents uncertainties in the luminosity. 
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luminosities lie between the hybrid/ AGN F03 ULIRG sam- 
ple and NGC6240 (Vignati et al. 1999), in which the X-ray 
emission is associated with a binary AGN system. The hard 
X-ray luminosities of the remaining sources are consistent 
both with the RCS03 and F03 starburst galaxies, as well 
as with the starburst-dominated local ULIRG Arp220 (Iwa- 
sawa et al. 2001). 



4 OPTICAL SPECTRA 

Figure 7 shows characteristic spectra for objects in the 
SB, LIRG and ULIRG luminosity classes (<z>~0.2, 0.5 
and 0.8, respectively). The DEEP 2 spectral range, ~6500- 
9000 A, permits the following principal lines to be visible: 
[NeIII](A3869), [Oil] (A3727) and &5(A4102) at 0.6 < z< 1.8, 
Ha (A6563) and [Nil] (A6584) at z < 0.3, H/3(A4861) and the 
[OIII] doublet ( A A4959, 5007) at 0.3<z<0.8. Emission line 
measurements (table 4) were made by using single Gaus- 
sian fits, after fitting the local stellar continuum with a 
linear combination of early-type and A-star spectral tem- 
plates, to approximate the effects of absorption and contin- 
uum shape on the line strengths. This is particularly im- 
portant for faint lines, such as [Nelll], which lie in spectral 
regions with strong continuum sub-structure. H<5 equivalent 
widths (EWs) were derived using the model-fitting method 
of Harker et al. (in prep.): for the Balmer lines, the emis- 
sion component, which comes from HII regions, follows a 
well-defined decrement dependent only on the extinction in 
the ionized gas. The EWs of the absorption line compo- 
nent, determined by the mean stellar populations, are also 
well-behaved across a range of star formation histories. By 
simultaneously measuring H/3 and H<5 in a spectrum, we are 
able to self-consistently estimate both emission and absorp- 
tion line EWs for the z ~ 0.6-0.8 subset of 70 fim and control 
galaxies. In a similar way, H/3 emission lines were corrected 
for the underlying H/3 absorption component. 

The diagnostic properties of emission lines, particularly 
useful when expressed in terms of line ratios, have found 
extensive use in determining conditions in the ISM, stellar 
ages, metallicity and the strength of the ionisation field. As 
means of separating nuclear and stellar excitation processes, 
one of the most widely-used techniques, proposed by Bald- 
win, Phillips & Terlevich (BPT, 1981), relies on identifying 
sources with respect to their [OIII] /H/3 and [Nil] /Ha ra- 
tios. As our sample spans a large redshift range, no emission 
lines are common to the entire sample and as a result wc 
cannot take advantage of such a diagnostic; there are no 
spectra with all 4 lines simultaneously visible. Ha is promi- 
nent in the < z > ~ 0.2 SB population, also seen in the low 
redshift LIRGs, but redshifted out of the spectral range at 
z > 0.35, where it gets replaced by emission from H/3. The 
strong Ha and H/3 signatures in the galaxies' spectra are in- 
dicative of rapid ongoing star-formation as they are directly 
related to the flux of massive stars ionising the ISM (e.g. 
Kewley, Geller & Jansen 2004). We find no indication of 
Doppler broadened (1000 < FWHM < 10000 kms" 1 ) hydro- 
gen lines suggestive of the presence of the broad line region 
(BLR) around an active nucleus. All sources with Ha in the 
spectral range also show the [Nil] feature, however the latter 
is significantly weaker (e.g. see Fig. 7). We also identify the 
4000 A spectral break (D4000) and the higher order Balmer 



absorption line H5, associated with the age of stellar popu- 
lations, in particular the presence of BAF class stars. D4000 
was defined by Bruzual (1983) as the ratio of the average flux 
density in the 4050-4250 A and 3750-3950 A band, with the 
narrower equivalent (D„4000: 4000-4100 A to 3850-3950 A) 
defined by Balogh et al. (1999). The existence of a break 
manifests a change in stellar opacity due to the accumulation 
of a large number of absorption lines from multiply ionised 
metals in the atmospheres of hot, metal-rich stars. Here we 
calculate D n 4000 as it is less sensitive to extinction, hence 
more appropriate for our sample and investigate its strength 
in relation to the Balmer absorption line H<5, both visible 
in the ~ 0.7-1.2 redshift range. However, as calculating H5 
EWs requires reliable H/3 measurements, we only have H<5 
absorption line information for the 8 objects in the sample 
for which both lines are in the accessible range (6 LIRGs and 
2 ULIRGs). Moreover, we are able to determine D n 4000 for 
only 7 out of 22 sources. For the remaining objects, we were 
not able to retrieve D n 4000 measurements due to one or 
a combination of the following: one of the D„ 4000 bands 
being partially blanked by bad pixels, emission line contam- 
ination from [Nelll], there being no evidence of a break or 
lack of the reliable spectrograph throughput corrections. 

Various common forbidden lines can also be seen in the 
spectra of the 70 /im sample, originating either in interstel- 
lar low-density PhotoDissociation Regions (PDRs) or the 
narrow line region (NLR) in the vicinity of an AGN. The 
luminosity and ionisation state of these lines is related to 
the strength of ionising UV flux. More specifically, high ion- 
isation states such as [NeV], which we detect in one object, 
are a consequence of a hard radiation field and hence unam- 
biguous indicators of the presence of an AGN. The strong 
[OIII] doublet appears in the 0.3< z <0.8 spectral range. The 
70 fim sources display various strengths of the [OIII] , [Oil] , 
[Nelll] and H /3 features and for the purpose of our study we 
are able to examine the [OIII] /H/3 and [NeIII]/[OII] ratios 
for 33 70 fim sources in total. 

In sections 4.1, 4.2 and 4.3 we examine the [OIII]/H/3 
and [NeIII]/[OII] ratios, as well as the behaviour of the H<5 
absorption line and its relation to the 4000A break and com- 
pare with equivalent measurements for the control galaxies 
matched to each 70 /xm source. Note, that for all calculations 
we rely on flux ratios between features at similar wavelength, 
so the need for extinction correction is eliminated. 

4.1 The [OIII] /H/3 ratio 

In traditional classification schemes (e.g. Ho, Filippenko & 
Sargent 1997) the [OIII]/H/3 ratio has been used to dis- 
tinguish starbursts and Seyferts (e.g. Baldwin, Wampler & 
Burbridge 1981; Veilleux & Osterbrock 1987), with higher 
ratios found in hybrid AGN/Starburst galaxies (e.g. Caputi 
et al. 2008) and a typical [OIII] /H/3 > 3 seen in Seyfert galax- 
ies. Here, we are able to measure [OIII] /H/3 for 20 LIRGs 
and 8 ULIRGs from the 70 /im sample and associated con- 
trol sources (Fig. 8). For the remaining sources one or both 
of these lines are out of range. 

The LIRGs appear to have approximately the same 
log [OIII] /H/3 mean value as the control galaxies (Fig. 8, 
top panel), roughly at -0.25, consistent with the line emis- 
sion being linked to stellar processes, apart from 2 ob- 
jects with [OIII]/H/3>5. Five LIRGs have the H/3 emis- 
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Figure 7. Characteristic spectra for the 70 fim sample, 3 SBs, 3 LIRGs and 3 ULIRGs. 



sion line out of range, but the [OIII] marginally in range, 
however the latter is weak and not measurable in all but 
one object. In terms of the control sample, a tail-off to 
high [OIII]/H/3 ([OIII]/H/3>3) is evident, possibly due to 
Seyfert-types which are likely a few percent of all galaxies at 
these redshifts. With respect to the ULIRGs (Fig. 8, lower 
panel), there is more spread in the two distributions, but 
with the 70 sample having higher mean [OIII]/H/3. Half 
the ULIRGs have [OIII]/H/3 < 3, below the rough dividing 
line between AGN and star-forming galaxies and half have 
3< [OIII]/H/3 < 10 typically seen in AGN. The 5 LIRGs and 
ULIRGs with [OIII]/H/3>3 ratios are considered in subse- 
quent work on the AGN fraction (section 5). 

4.2 The [NeIII]/[OII] ratio 

The [NeIII]/[OII] ratio is a metallicity indicator, where the 
negative correlation with metallicity arises due to the atmo- 
spheres of more metal-rich stars absorbing a greater frac- 
tion of ionising photons, softening the radiation field and 
decreasing the abundance of the higher ionisation states (e.g. 
Nagao, Maiolino & Marconi 2006). [NeIII]/[OII] is also an 
ionisation discriminator with AGN showing higher ratios on 
average, typically ~ 0.4 for type 2 and >1 for typel (e.g. 
Nagao et al. 2002). 

The [NeIII]/[OII] ratio for our sample can only be exam- 
ined for the z>0.7 sources, so all, apart from two, of the 10 
sources with reliable measurements are ULIRGs (table 4 and 
Fig. 9). The majority of the control galaxies are assembled 



between -1.5<log [NeIII]/[OII]<-0.4, which is the range one 
would expect for distant star-forming galaxies (e.g. Perez- 
Montero et al. 2009), with 2 galaxies at log [NeIII]/[OII]>- 
0.5, which are possibly Seyferts or metal-poor. The 70 fim 
sources have -l<log [NeIII]/[OII]<-0.2 in two groups which 
peak at ~-0.7 and ~-0.3. The first group of 5 sources is in 
agreement with the average [NeIII]/[OII] ratio of the control 
sample. On the other hand, the peak at log [NeIII]/[OII]> - 
0.4, occupied by the remaining 5 70 /im sources, is degen- 
erate as it is consistent with both a type 2 AGN regime 
(e.g. Nagao et al. 2002) and a low metallicity regime (e.g. 
Perez- Montero et al. 2007). 

4.3 The H<5 absorption line and 4000 A break 

Figure 10 compares the H 8 absorption line for the 7 sources 
from the 70 fim sample which are in the right redshift range 
to have good line measurements, to their respective con- 
trol galaxies. The presence of higher order Balmer absorp- 
tion lines are indicative of recent star-forming activity (e.g. 
Dressier & Gunn 1983) , with the peak in H 8 absorption oc- 
curing at about ~ 0.5 Gyr, followed by a decrease during the 
subsequent 1-10 Gyr epoch. The sources in our sample have 
higher H 8 EWs than the bulk of the control sample, im- 
plying a younger mean stellar population. Nevertheless, H5 
EWs in the 5-10 range are age degenerate, corresponding to 
systems that are either very young (~ 10 Myr) or ~1 Gyr old 
(e.g. Kauffmann et al. 2003, hereafter K03). To settle this 
degeneracy we measure the narrow 4000 A break (D n 4000) 
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Figure 8. Distribution of the [OIIIJ/H/3 ratio for objects in the 
70 (im sample (solid lines) — top panel for the LIRGs and lower 
panel for the ULIRGs — and their corresponding control galaxies 
(black dot-dashed lines). 




Figure 9. Distribution of the [NeIII]/[OII] ratio for objects in the 
70 fim sample (red solid line) and their respective control galaxies 
(black dot-dashed line). 



Figure 10. The distribution in H<5 absorption line EW. Red solid 
lines: the 7 sources from the 70fim sample which are in the right 
redshift range for reliable measurements (see section 4). Black 
dot-dashed lines: their corresponding control galaxies. 
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Figure 11. Plot of H<5 EW vs D„4000. Red filled circles: the 
control sample, black diamonds: the 70 fim sample. 



break and examine its strength against the H<5 absorption 
EW in figure 11. The 4000 A break is a monotonic func- 
tion of age, small for young stellar populations and large for 
old, metal-rich galaxies. We find that both the control and 
70 fim sample lie where expected for star- forming galaxies 
with EW(H(5) > 1 and D4000 > 1, however their spread over 
the D n 4000-H<5 plane is substantially different: the control 
sample spans a large part of the plane, whereas the 70 )im 
sources have low D n 4000 (< 1.3) and most have high EW 
(H<5) (> 5) indicative of young stellar populations and con- 
sistent with a starburst mode of star formation (e.g. K03). 
Similar values of D„4000 and H 8 are also reported by Mar- 
cillac et al. (2006) and Caputi et al. (2008)) who examine 
distant IR-selected galaxies. Caputi et al. (2008) find them 
to have the smallest D n 4000 values amongst a large sample 
of sources from the COSMOlogical evolution Survey (COS- 
MOS), on par with the 70 fim sample occupying the lower 
end of the D n 4000 distribution of DEEP 2 control galaxies. 
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Figure 12. DEEP 2 spectra for the 8 AGN candidates, whose SEDs are shown in figure 13. 



5 THE AGN CONTENT OF THE 70 fim 
POPULATION 

5.1 The AGN fraction 

In order to calculate the AGN fraction, we use four diagnos- 
tics — hard X-ray emission, high (>3) [OIIIJ/H/3 ratio, the 
presence of [NeV] lines and a power-law near/mid-IR con- 
tinuum — with the requirement that at least one is satisfied 
for an object to be classed as an AGN candidate; see table 
2 and figures 12 and 13. 

The X-ray fraction in the 70 jj,m sample amounts to 
16 (±6) per cent: two objects are detected in the full X- 



ray band only (objects 40 and 84), whereas 8 have detec- 
tions in at least one other band (objects 43, 56, 59, 67, 
83, 93, 122 and 147). However, as star formation can also 
be a source of X-ray emission (see section 3) we only con- 
sider sources with hard X-ray detections to be AGN candi- 
dates — objects 59, 67, 83 and 122 — especially since for 
the remaining sources X-ray emission is consistent with the 
star-forming relations in figure 6. These 4 sources also have 
log L,2_iokcV > 42, which, as previously mentioned, usually 
indicates that hard X-ray emission is AGN-dominated. 

A power-law type near/mid-IR (3 < A < 20 /im) contin- 
uum and hence red colours in that part of the SED, are a 
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Figure 13. The X-ray to IR SEDs for the 8 objects in the sample which show at least one AGN signature; black filled circles: optical/IR 
data, black solid line: SK07 IR model SED templates (normalised to 24 /im); dotted line: radio quiet QSO SED from Elvis et al. (1994) 
(normalised to the absorption-corrected hard-band X-ray flux), dash-dot line: SED of AGN NGC7213 from Kuraszkiewicz et al. (2001) 
(normalised to the absorption-corrected hard-band X-ray flux), diamonds: hard and soft fluxes corrected for absorption, triangles: hard 
and soft fluxes not corrected for absorption (red for measured flux, blue for upper limit). When there is no X-ray detection the AGN 
SEDs are normalised to the hard X-ray upper limits. For object 165, which does not have an X-ray detection, but has a power-law 
near-IR continuum, the templates are normalised to both the hard X-ray upper limit and the near-IR continuum. 



consequence of emission from dust heated to near sublima- 
tion temperatures. As the strength of the stellar radiation 
field is not sufficient to cause this, the only stellar-related 
mechanism that could be responsible for a near/mid-IR ex- 
cess is dust stochastically heated by shocks from outflows 
and supernovae (e.g. Ho et al. 1989; Davies, Burston & Ward 
2002). However, studies have shown that the dust grains re- 



sponsible are of very small size and hence lifetimes of the or- 
der of only a few kyr (Dwek 1986), making this an unlikely 
explanation for a near-IR excess in a broad-band SED. A 
more likely scenario is that of AGN dust heating and hence 
direct emission from the torus (Neugebauer 1979; Elvis et 
al. 1994; Lutz et al. 1998; Sturm et al. 2000; Klaas et al. 
2001; Alonso-Herrero et al. 2006). There is one object in 
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Figure 14. The X-ray to IR SEDs of 6 objects with no AGN signatures (the LIRGs and ULIRGs in Fig. 7) and AGN-candidate object 
165; black upper limits: X-ray data; black filled circles: optical/IR data; black solid line: SK07 IR model SED templates (normalised to 
24 fim). The red open diamonds represent the panchromatic SED of NGC 1068 — photometry for the SED from Sandage (1973); Spinoglio 
et al. (1995); Stickel et al. (2004); Cappi et al. (2006) and Howell et al. (2007). The green filled squares represent the panchromatic SED 
of NGC 6240 — photometry for the SED from Allen (1976); Soifer et al. (1989); Iwasawa & Comastri (1998); Spinoglio et al. (1995); 
Klaas et al. (2001) and Lutz et al. (2004). The SEDs of NGC 1068 and NGC 6240 are normalised to the SK07 template at the IRAS 
60/im photometry, in order to be consistent with our original sample selection at 70^tm. Note that for both local galaxies including the 
soft X-ray flux is not relevant, as it is dominated by emission lines from photoionized plasma, with the underlying X-ray continuum 
completely absorbed. The hard X-ray flux that we plot here is observed and not corrected for absorption. 



the 70 /jm sample with a power-law near/mid-IR continuum 
and red colours in all IRAC bands (3.6-8 (im) — object 165 
(see SED in Fig. 13). It also has a high (>3) [OIII]/H/3 
ratio; see table 2. The remaining AGN candidate sources 
have starburst-type SEDs, although this does not exclude 



an AGN contribution to the near/mid-IR. In terms of AGN 
spectroscopic signatures, object 67 has an evident [NeV] line 
and objects 84, 114, 141, 122, 165 have high [OIII]/H/3 ratios 
(objects 67 and 122 are also hard X-ray emitters). Object 83, 
one of the hard X-ray emitters, has a spectrum that is almost 
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devoid of emission lines, which implies complete obscuration 
or absence of the narrow line region and could be an exam- 
ple of a host-obscured source (e.g. Martinez-Sansigre 2006), 
or a source with near 4n torus covering factor (Imanishi et 
al. 2007). 

In total, 8 out of 61 sources in the 70 pm sample show 
evidence for an active nucleus by satisfying at least one di- 
agnostic, making the AGN fraction 13 per cent, or 11 and 
23 per cent when the LIRGs and ULIRGs are considered 
seperately. Note that this is the fraction of sources which 
host an AGN, not the fraction of sources powered by AGN 
and it is likely a lower limit, as it is possible that a number of 
AGN without spectroscopic or photometric signatures could 
be missed — see section 5.3 for a discussion. 

5.2 The AGN contribution to the energy budget 

For the 8 sources which are found to host an active nucleus, 
we investigate the AGN contribution to the total energy 
budget with the aid of SED templates; figure 13 shows the 
sources' panchromatic SEDs with X-ray fluxes in the hard 
and soft bands, B, R, I fluxes from DEEP 2, J and K fluxes 
where available, IRAC and MIPS data and the SK07 IR 
template best matched to each object (see section 2.1). We 
use Webpimms, a NASA High Energy Astrophysics Science 
Archive Research Center (HEASARC) tool, to estimate ab- 
sorbing column densities (Njj) and correct the fluxes for 
absorption, via our estimated hardness ratios and a r=1.9 
photon index. These are shown in table 4 where the upper 
and lower values for the column density are calculated using 
the upper and lower values of the hardness ratios. We use the 
average SED for a radio quiet QSO from Elvis et al. (1994) 
and that of a lower luminosity, lower Eddington-ratio AGN, 
NGC 7213 from Kuraszkiewicz et al. (2003) which lacks the 
'big blue bump'. Assuming that hard X-rays are a measure 
of AGN power, both templates are normalised to the absorp- 
tion corrected hard X-ray flux or the hard-band upper limit 
if the object either has no X-ray detection or it is detected 
only in the full-band. In order to separate the AGN and host 
galaxy SED components we use templates which have little 
or no contribution from the host galaxy and therefore rep- 
resent integrated energy from the AGN. Note that because 
these templates represent unabsorbed AGN, they would not 
be expected to match our observed spectrum in the opti- 
cal/UV according to the absorbing columns we calculate. 
Nevertheless and assuming the unified scheme holds (Urry & 
Padovani 1995), the difference between an unabsorbed and 
an absorbed AGN SED should be far more obvious in the 
soft X-ray/EUV/UV and optical, whereas the contribution 
from the torus in the mid and far-IR, where it is optically 
thin, should vary much less with orientation. Consequenctly, 
since our aim is to calculate the integrated AGN energy by 
normalising the templates to our sources' photometry, it is 
not appropriate to use absorbed AGN templates where a 
large part of the spectrum would be partially or completely 
wiped out. This is particularly important for higher lumi- 
nosity AGN, where most energy comes out in the UV and 
optical. Moreover, unabsorbed templates allow us to evalu- 
ate the strength of the stellar bump in comparison to the 
AGN in the 1-2 /im region. 

The SEDs of all 8 sources (Fig. 13) are characterised by 
a far-IR bump significantly more luminous than the AGN 



templates in that part of the spectrum. For objects 84, 114 
and 141 the entire optical/IR SED is displaced from the 
AGN templates and is 2-4 orders of magnitude more lumi- 
nous, an indication that the AGN is of low luminosity and 
therefore unlikely to contribute substantially to the galaxies' 
energy budget. For the hard X-ray emitters, objects 59, 67, 
83 and 122, the AGN templates agree better with the mid-IR 
SEDs, indicating that some part of the galaxies' near/mid- 
IR luminosity must be coming from the AGN. Object 165 
is not detected in the X-rays, but it is classed as hosting an 
AGN from spectral signatures and mid-IR continuum slope. 
Accordingly it is more appropriate to normalise the AGN 
template to the near/mid IR continuum, as the power-law 
SED slope in that region is characteristically emission from 
the torus, showing that even if the AGN is not detected 
in the X-rays, it could be contributing substantially to the 
source's energy budget. 

We examine the possibility that the far-IR bump re- 
sults from the high frequency light from the AGN being 
absorbed and re-radiated into the IR. Assuming that X-ray 
emission scales with AGN power, then to power a LIRG or 
ULIRG (Lm : = 10 n -10 13 L ) and hence dominate emission 
in the far-IR, we would expect an AGN luminosity of at 
least a few times 10 44 erg/s and assuming that a few per 
cent is emerging in the X-rays (e.g. Elvis et al. 1994), an 
X-ray luminosity of at least 10 43 erg/s for the low luminos- 
ity IR sources and at least 10 44 erg/s for the more luminous 
LIRGs and ULIRGs. None of our objects satisfy these crite- 
ria. In fact, quite simply integrating under the SEDs, reveals 
that in all cases the IR bump (8-1000 fim) has between 1-3 
orders of magnitude more energy than the AGN/QSO tem- 
plates (defined as the Lir/Lagn ratio). It is possible that 
the column densities, and hence absorption corrections, de- 
rived from our hardness ratios could be underestimated if 
there is a contribution from scattered radiation in the soft 
X-ray band for heavily absorbed sources. In order to ad- 
dress this, we estimate a maximum value for Lir/Lagn by 
normalising the AGN template to our 8 /jm IRAC photom- 
etry. We find that the IR SED overtakes emission from the 
AGN by factors of 6-90 and hence conclude that a) with 
the exception of object 165, where Lir/Lagn is less than 
an order of magnitude, for the remaining objects the AGN 
contributes less than 10 per cent to the IR energy budget 
and b) all sources are primarily powered by star-formation 
— see table 2. 

5.3 Obscured AGN 

Our results indicate that low luminosity 
(L2-iofceV r ~ 10 42 erg/s) AGN with moderate obscura- 
tion (Nh ~ 10 23 cm~ 2 ) can be detected in our survey up to 
z~0.5. Active nuclei obscured by higher column densities 
are potentially missed even if more luminous, implying 
that the AGN fraction that we calculate is a lower limit. 
Nevertheless, the hardness ratios we calculate for our 
stacked sources (table 1) are soft (HR<-0.2) indicating a 
minor (if any) contribution from obscured AGN to the 
X-ray signal, making our AGN fraction relatively robust 
- this is further confirmed by the stacking signal from 
the same sample of sources minus the 8 AGN candidates 
(see table 3). In addition, even without X-ray detections, 
it is unlikely that we would mis-classify AGN-powered 
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Table 2. The 8 objects identified to host AGN. Columns 1-6: object ID, infrared luminosity class, SED type, hard X-ray luminosity, the 
[OIIIJ/H/3 ratio, a visible NeV line. Lir/Lagn is the ratio of the integrated energy in the 8-1000 (im spectral region to the integrated 
energy in the QSO template from Elvis et al. (1994) when normalised to the absorption-corrected hard X-rays or X-ray upper limit 
(column 7) and the IRAC 8 /mi photometry (column 8). The latter represents the maximum AGN contribution as it is constrained by 
our infrared photometry. Column 9 identifies the primary energy source of the galaxy. Given the power-law-type SED of object 165, 
the AGN likely has a substantial contribution to the galaxy's energy budget, however the ^ir/^agn ratio shows that the starburst 
component is still the primary energy source. 



ID 


Lift class 


SED type 


log L 2 _iokev(erg/s) 


[OIII]/H/3 


[NeV] 


Lia/Lagjv (X) 


Lir/Lagn (8) 


Powered by 


59 


LIRG 


starburst 


42.9 


7 


out of range 


23 


10 


starburst 


67 


ULIRG 


starburst 


43.3 


out of range 


yes 


59 


40 


starburst 


83 


LIRG 


starburst 


42.6 


too weak 


out of range 


42 


13 


starburst 


84 


ULIRG 


starburst 


not detected 


9 


out of range 


2448 


90 


starburst 


114 


ULIRG 


starburst 


not detected 


3 


out of range 


640 


27 


starburst 


122 


LIRG 


starburst 


43.4 


out of range 


out of range 


14 


13 


starburst 


141 


LIRG 


starburst 


not detected 


6 


out of range 


192 


26 


starburst 


165 


ULIRG 


power-law 


not detected 


8 


out of range 


550 


6 


starburst /AGN 



Table 3. Table of stacked fluxes (erg/s/cm 2 ) and hardness ratios for the 70 (im sample minus the 
8 AGN candidates (see section 5.1). They are split into 3 luminosity classes (starbursts, LIRGs 
and ULIRGs; see section 2.1 for details on the sources' infrared properties). The second column 
indicates how many objects were used in the stacking — some were removed because they were 
too close to a source which could have contaminated the signal or they were associated with an 
X-ray source. The SB group does not have a formal detection in the hard band, so we quote the 
3cr upper limit. The SB hardness ratio is derived using the counts in the soft band and 3cr upper 
limit counts in the hard band. 

Class Objects used fo.5-2 h-io HR 

Starbursts 8 out of 9 8.77E-17 (8er) 2.35E-16 (3a upper limit) -0.3 (upper limit) 
LIRGs 26 out of 31 5.93E-17 (lOer) 1.6c-16 (4er) -0.3 

ULIRGs 9 out of 13 9.82E-17 (10cr) 2.33E-16 (3<r) -0.4 



70 pirn sources as starburst-dominated. As outlined in 
sections 5.1 and 5.2, AGN and starburst SEDs are highly 
divergent in the near/mid IR and therefore if one of our 
L>1O 11 L sources were to host a L2-iofc e v~ > 10 44 erg/s, 
haevily obscured AGN which could potentially power the 
galaxy but is not detected in the X-rays, its SED would 
show a strong near/mid-IR continuum: hot dust emission 
from the torus. 

To demonstrate this further we compare the SEDs of 6 
objects with no AGN signatures (the LIRGs and ULIRGs 
in Fig. 7) and AGN-candidate object 165 to the SEDs of 
the local IR galaxies NGC 6240 and NGC 1068 known to 
host compton thick AGN (Fig. 14). The diverse nature of 
NGC 6240 and NGC 1068 is unambiguously mirrored in their 
panchromatic SEDs, which are in agreement over the optical 
and far-IR parts but diverge significantly in the near/mid- 
IR. NGC 6240 is a SB/ AGN hybrid with a binary compton 
thick (~10 24 cm -2 ), ~10 44 erg/s QSO, but the galaxy is nev- 
ertheless primarily powered by a starburst (e.g. Vignati et al. 
1999; Komossa & Schultz 1999; Klaas et al. 2001), whereas 
NGC 1068 is a type-II Seyfert (L 2 _iofe e v >10 44 erg/s) with 
an additional far-IR starburst component (e.g. Le Floc'h 
2001; Spinoglio et al. 2005). We normalise their SEDs to the 
SK07 template at the IRAS 60/im photometry, in order to 
be consistent with our original sample selection at 70pim. 
We do not consider the soft X-ray flux, which has limited 
diagnostic value for such heavily obscured AGN embedded 
in star-forming galaxies. The hard X-ray flux we show in the 
SEDs of NGC6240 and NGC1068 is observed and hence not 
corrected for absorption. 



We see from Fig. 14 that if there are any sources in 
the 70 /im sample with IR to X-ray (L60/jm/L 2 _iofcev) ratios 
similar to NGC 6240 and NGC 1068 then they would on av- 
erage be detected in our X-ray survey. However, even if such 
compton thick AGN are missed, the near/mid-IR part of the 
SED would show unambiguous signatures of their presence 
if they were to be significant contributors to the galaxy's en- 
ergy budget. Not surpsisingly, apart from object 165, whose 
SED is better represented by that of NGC 1068, the remain- 
ing 70 fim sources are in agreement with the starburst-type 
SED of NGC 6240. 



Table 4: The 70 /im sample (61 objects): IDs from the 70 /jm EGS survey 
and the DEEP 2 survey, redshift, R-band magnitude (AB), total infrared 
luminosity, X-ray flux in the soft and hard bands (xl0~ 16 ergs -1 cm 2 ), 
hardness ratios, column density (cm -2 ) with the upper and lower limits, 
line ([Oil], [Nelll] H/3, [OIII]) fluxes (log W/m 2 ), H5 absorption line 
equivalent width (A) and the reduced 4000 A break. '(1)' indicates a 
3a upper limit for the fluxes. Column densities were calculated with 
r=1.9. '(1)' in the N_h column indicates that with a T=1.9 the estimated 
column density is lower than the Galactic Nh of ~ 2 x 10 20 . The upper 
and lower limits of Nh are calculated from the upper and lower values 
of the hardness ratio. 



ID 


DEEP ID 


z 


R mag 


Lir x 1O 1U (L ) 


fo.5-2 


f2-10 


HR 




{[Oil] 


f[NeIII] 


fff/3 


f[o///] 




n 4000 


STARBURSTS 






20 


14017095 


0.33 


19.70 


8.72±0.11 






















46 


13056925 


0.22 


19.72 


9.24±0.04 






















51 


13058203 


0.22 


18.59 


4.37±0.23 






















56 


13041622 


0.20 


18.70 


9.84±0.49 


9.04±i;I 5 


13.7(1) 


-0.48±8i 7 
















81 


13009920 


0.23 


19.65 


5.38±0.38 




















118 


12008271 


0.24 


20.30 


3.77±0.30 






















125 


11051641 


0.25 


20.72 


8.28±0.12 






















144 


11033888 


0.19 


19.35 


3.62±0.18 


















_ 


_ 


155 


11020745 


0.17 


18.60 


2.74±0.30 


















_ 


_ 


LIRGs 






33 


14012882 


0.34 


19.80 


18.72±2.87 
















-19.23 






37 


13063596 


0.36 


19.49 


45.23±0.46 














-18.86 


-19.43 






40 


13063597 


0.30 


18.36 


41.96±2.10 


1.75(1) 


10.4(1) 


-o.24±8:i 


5.4X10 21 J,, 6 * 1022 








-19.70 






41 


14007624 


0.45 


20.88 


80.59±7.03 














-19.61 


-20.06 




- 


43 


13063920 


0.56 


20.73 


32.42±3.09 


i.i3±8:§ 3 7 


2.23(1) 


-o.66±g : i 


/i\ 4xl0 21 
W (1) 






-19.03 


-19.39 






55 


13049741 


0.67 


21.56 


84.70±5.68 














-19.33 


-20.15 


5.52 


1.19 


59 


13050479 


0.69 


22.33 


82.16±2.77 


4.99±i;i 


106±J 2 ; 7 


o.6i±8:i ,s 


i i v m 23 i-isxio 23 

1.1 X 1U 8xl0 22 




-19.86 


-19.70 


-18.88 


0.25 


1.29 


68 


13034644 


0.35 


19.65 


33.96±1.50 














-19.12 






72 


13035302 


0.39 


20.44 


13.77±1.18 
















-20.07 






76 


13017944 


0.73 


22.62 


30.96±2.16 














-19.54 


-19.68 


4.45 


1.14 


78 


13026857 


0.37 


20.27 


35.92±0.15 














-18.91 


-19.58 






83 


13019240 


0.75 


22.55 


43.08±11.08 


5.53±i:ii 


35.5±7;g9 


o.i5±8:il 


4 2 x 10 22 B - 7xl ° 22 

4.Z X 1U 2 6xl0 22 


-19.99 


-20.73 




-19.74 






85 


13003379 


0.43 


19.92 


25.80±0.92 












-18.82 


-19.54 






94 


12027947 


0.55 


21.69 


76.92±3.03 














-19.17 


-19.77 






95 


12027969 


0.73 


22.01 


90.08±29.67 






















98 


12023868 


0.29 


20.10 


10.82±1.19 






















100 


12024073 


0.30 


19.33 


23.56±1.57 






















103 


12029058 


0.25 


19.02 


10.90±0.15 






















110 


12015978 


0.42 


21.29 


33.28±2.03 














-20.42 


-20.42 







3 

as 
3 

2 



03 
O 

s 

CO 



03 

03 
C3 

ro 

Ei. 

a 
^ . 

03 
03 

CO 

S3 
-i 
o- 

S 

2 

o 

-5 
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ID 


DtLtLr 1JJ 


z 


R mag 


Lin x 1U (Lq ) 


£ 

10.5-2 


r 

12-10 


T IT) 




£ 

Ho 1 1] 


r 

L[NeIII] 


c 

Iff ,3 


r 


hj\\ HS 


Ti a nnn 
JJ4UUU 


1 in 

112 


LZUZU til 


n o a 


i n oo 


07 £?0_1_7 70 

2/.bo± ( . to 
















m AG 

-iy.4b 






llo 


i onnQ070 
12UU02 ( 6 


n oo 
U.oo 


on no 
2U.U8 


OK 1 O-l-O 1 7 

2o.l2±2.1 / 














1 Q 77 


i n O/i 
-iy.24 






i on 
12U 


i oni oi c7 
12UlolD( 


n c7 


oi a 1 
21.61 


on i o_i_o 7^ 
o0.1o±o. to 














i n oo 

-iy.oo 


i n i o 

-iy.12 






121 


i onnom A 
12U08914 


n to 
U. (0 


oo n7 
22. U 1 


c 7 no_i_o en 
o/.U8±8.oO 








' X 1U 6xl<) 22 

- 


i n on 
-19.29 


i n c a 

-iy.b4 


i n i 7 

-iy.i t 


i n /i n 

-iy.4y 


7 77 


1.1 


122 


12004450 


0.28 


19.42 


11.80±2.70 


3.56±l:i 


„ _~ ,17 

172±iO 


0.78±g;8f 


- 


- 


- 


- 




- 


126 


11051657 


0.42 


21.16 


27.97±2.23 


- 


- 


- 


- 


- 


- 


- 




- 


135 


11046414 


0.48 


20.92 


57.20±0.28 


- 


- 


- 




- 


- 


- 


- 




- 


137 


11046419 


0.48 


20.07 


67.64±0.39 


- 


- 


- 


- 


- 




-18.61 


-18.75 




- 


1 oo 

loo 


110^9469 


U.42 


OI 1 o 

21. lo 


i o c /i 07 
12.o4±b.2 I 














i n oo 

-iy.2o 


inn 

-iy.y 






1 /I n 

14U 


1 1 A/l 7nQ£? 

lllKWUoo 


n /i r 
U.4D 


i n /i o 
ly.42 


07 no_i_ £ fio 
2/.Uo±0.02 






















141 


i i i m i 77 
111U11 f f 


n 7n 
U. 1 U 


OO 7/1 

22. ( 4 


0£ 7£_1_1 O/l 
00. / 0±1.04 












-2u.uy 


i n £Q 

-ly.oo 


1 Q 7n 

-lo. t\> 


1 A 1 A 

1U.14 




146 


11033882 


0.46 


20.05 


19.64±2.18 








/i\ 1.8X10 22 
W (1) 














147 


11026910 


0.67 


22.31 


85.49±14.92 


i.53±8:S 


2.58(1) 


-o.5±8:l 9 






-19.74 


-20.61 






149 


11027507 


0.51 


20.34 


46.83±7.07 














-18.52 


-18.45 






150 


11034350 


0.74 


21.97 


48.99±8.80 














-19.44 


-19.82 


0.O4 


1.1 / 


152 


11027541 


0.37 


19.98 


54.70±0.37 














-19.19 


-19.42 






ULIRGs 






24 


14011285 


0.76 


22.07 


153.57±23.59 










-19.82 




-20.28 


-19.89 






26 


14017213 


1.29 


23.68 


796.20±51.18 










-19.38 


-20.23 










58 


13041891 


0.69 


22.29 


196.78±42.96 








Q 7 ^ 1 n 23 3.9X10 23 
6 -' X 1U 2.8X10 23 






-19.67 


-20.35 


^ Q1 


1 Oft 
1 .Zo 


67 


13034619 


0.96 


22.11 


235.69±2.43 


0.3(1) 


55.2±xcl! 


o.93±8:S 


-19.17 


-19.56 










70 


13026142 


0.78 


21.30 


223.68±56.86 








-19.62 


-20.32 










77 


13036021 


0.42 


19.16 


101.04±0.60 














-18.96 


-19.52 






82 


13027191 


1.04 


23.36 


275.11±44.43 








6 3 x 10 22 1 ■ lxll)23 

0.0 X 1U 2 3xl0 22 


-19.34 


-19.73 










84 


13019950 


0.77 


21.97 


307.45±64.62 


0.92(1) 


4.2(1) 


o.32±8:f 


-19.20 


-19.50 


-19.34 


-18.39 






88 


13019982 


0.78 


21.51 


254.05±10.51 








-19.16 


-19.51 


-19.78 


-19.40 






93 


13004291 


1.20 


22.57 


488.97±74.57 


i-59±8:l? 


3.82(1) 


-o.34±8:i 


1 X 10 2 ™ 22 


-19.54 


-20.50 






- 


- 


101 


12028577 


0.82 


23.63 


100.84±12.61 
















- 


- 


114 


12007915 


0.66 


20.83 


128.59±2.51 














-19.66 


-19.18 






124 


11051625 


0.85 


21.31 


173.78±43.18 










-19.2 








- 


- 


134 


11038951 


0.75 


21.87 


107.38±1.00 












-20.30 


-19.99 


-19.89 






158 


11020790 


0.93 


22.08 


795.07±163.23 










-19.58 












165 


11014633 


0.77 


21.71 


126.03±9.91 










-19.02 


-19.63 


-19.35 


-18.46 


5.76 


1.14 


172 


11015249 


0.81 


23.11 


228.37±1.09 
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6 SUMMARY AND CONCLUSIONS 

We have investigated the nature of 61 70 /mi-selected 
sources, by exploiting data from Spitzer, Chandra and the 
Keck telescopes, the latter as part of the DEEP 2 survey. In 
previous work we showed that the energy budget of 70 /im- 
selected galaxies is defined by emission in the infrared, with 
all z>0.1 galaxies having luminosities above Lin > 10 10 Lq. 
Our results revealed that 70 /im populations are predomi- 
nantly comprised of LIRGs at an average redshift of z~0.5 
and star formation rates of the order of ~ 100 Mq /jr. In 
addition, we found that the strong majority are identified 
with a starburst-type SED: a well-defined optical/near-IR 
stellar bump, followed by an inflection at 2< A rest < 6/tm 
and a sharp increase in flux at infrared wavelengths. The 
aim of the work described in this paper was to examine the 
X-ray and spectral properties of these sources and, by build- 
ing on previous results, develop a more detailed picture of 
the inherent nature of 70 /tm-selected galaxies. 

In order to place our sample in the context of other 
galaxy populations in the DEEP 2 photometric and spectro- 
scopic survey, we evaluated its spectral properties against 
a control sample assembled from DEEP 2 objects of simi- 
lar redshifts and optical colours (section 2.3). The clearly 
identifiable stellar bump in our sources' SEDs, places them 
in the optically bright (-18<M s <-23) and red (0.5 < U- 
B<1.5) regime in the colour-magnitude diagram (CMD). 
Our results showed that although there is agreement be- 
tween some control and 70 /im galaxies, there were also evi- 
dent differences both in terms of their AGN and stellar con- 
tent. The [OIII]/H/3 distributions of 70 /mi LIRGs and re- 
spective control galaxies were in agreement, with only a few 
outliers, implying that the optically-determined AGN frac- 
tions are broadly consistent in that redshift and U-B/Mb 
range. In contrast, there was a more evident shift towards 
higher [OIII]/H/3 for the ULIRGs than their respective con- 
trol galaxies. In terms of the [NeIII]/[OII] ratio half of the 
70 /jm sample was found to agree with the comparison sam- 
ple, whereas half was offset towards higher values. As dis- 
cussed in section 4.2, high [NeIII]/[OII] ratios can be repre- 
sentative of low metallicity systems, but also a harder ioni- 
sation field, albeit not unambiguously. In terms of H<5 EWs 
and the H<5-D n 4000 relation, the spread in values for the two 
samples show little overlap. The 70 /im population seems 
to be at a different evolutionary state to the control sam- 
ple, undergoing intense starburst episodes and characterised 
by younger stellar populations. This might not be surpising 
since, as mentioned earlier, the place of the 70 /mi sample on 
the CMD diagram could have been affected by extinction. 

Our main results are: 

• A total of 8 out of 61 sources (13 per cent of the sample) 
show at least one AGN signature and are therefore classed as 
hosting an active nucleus. The AGN fraction is 0, 11 and 23 
per cent for the SB, LIRG and ULIRG groups separately, ro- 
bust for L 2 -iofceV > 10 42 erg/s and N H < 4xl0 23 cm~ 2 , but 
potentially an overall lower limit due to a number of lower lu- 
minosity, more heavily obscured AGN that could be missed. 
[Note that by AGN fraction, we simply refer to the num- 
ber of sources which host an AGN]. The increase with total 
infrared luminosity is consistent with both local and high 
redshift studies (e.g. Lutz et al. 1998, Fadda et al. 2002, 
Franceschini et al. 2003). 



• 7 out of the 8 sources which show signatures of an active 
nucleus, have SEDs which are starburst-type with the inte- 
grated infrared emission ~ 1-3 orders of magnitude higher 
than the integrated AGN emission. One object has a power- 
law type SED, with a strong near/mid-IR continuum which 
implies that the AGN is energetically important. Obscured 
AGN which are possibly missed when estimating the AGN 
fraction, are not likely to play a role when estimating the 
fraction of AGN-dominated sources. We estimate that for 
an AGN to power a LIRG or ULIRG, it must have an X- 
ray luminosity of at least 10 43 erg/s and up to 10 45 erg/s. 
If the central black hole activity is powerful enough to con- 
tribute substantially or even dominate the galaxy's energy 
budget, unambiguous AGN signatures would emerge in the 
near/mid-IR part of the SED, even if missed in the X-rays. 
As a result, we conclude that all sources in the 70 /im sam- 
ple are primarily powered by star-formation — this includes 
object 165 which has an energetically important AGN, as 
its far-IR SED is more luminous than the AGN component. 

• The X-ray detection fraction of the control sample is 
4.2(±1.6) per cent, whereas for the 70 /mi sources it is 
16.4(±6.4) per cent. In terms of the latter, about 50 per 
cent comes from the objects we identified as AGN, showing 
a strong link between the X-ray and AGN fractions. This 
is to be expected, as the EGS Chandra survey is not deep 
enough to be sensitive to most of the star-forming sources 
in our sample. The higher X-ray, and hence, AGN incidence 
for the 70 /mi sample is at the ~2a" significance level. It pos- 
sibly relates to the higher estimated average stellar mass, 
as black hole mass is expected to scale roughly with galaxy 
mass, however, it could also hint at a relation between AGN 
activity and star-formation (e.g. see also Kauffmann et al. 
2003 and Silverman et al. 2009 find a number of AGN that 
reside in infrared-selected star- forming hosts). It remains to 
be seen whether the higher AGN fraction in 70/xm sources 
is due to the higher stellar mass or the more intense star- 
forming activity or both. 

• The AGN incidence of 13 per cent that we estimate is 
significantly lower than what has been previously observed. 
Previous studies have shown that, although star-formation 
is the dominant energy source in LIRGs and ULIRGs, a 
significant fraction of them host bolometrically important 
AGN (e.g., Solomon et al. 1997; Vcilleux, Sanders & Kim 
1997; Genzel et al. 1998; Downes & Solomon 1998; Scoville 
et al. 2000; Soifer et al. 2001; Fadda et al. 2002; Brand et 
al. 2006). Given this discrepancy, one would naturally ask 
whether there is a bias with respect to the wavelength of se- 
lection. The current status of far-infrared astronomy implies 
that the majority of studies on the AGN content of infrared 
galaxies have been based on populations selected in the mid- 
IR with the MIPS 24 /mi band. As discussed in section 5, 
there are fundamental differences between AGN and star- 
burst SEDs and although entirely AGN or starburst dom- 
inated objects are the extremes of this range, objects with 
AGN /starburst components of comparable contributions are 
rare and as a result sources will fall in one of the two cate- 
gories. Our results have shown that a far-IR bump and an 
SED that peaks longward of ~50 /im is strong evidence for 
the source being starburst-dominated, as AGN-dominated 
sources peak at much shorter (near/mid-IR) wavelengths 
and are weak far-IR emitters (see also Alonso-Herrero et al. 
2006). Consequently, one would expect a higher AGN inci- 
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dence in 24 /im populations, as a) it probes emission from 
warm/hot dust and b) at the mean redshift of the popu- 
lation (z~0.8) it corresponds to A re st~13/im, where there 
is little relevance to ongoing star formation. It seems that 
setting the selection at 70/im, at least down to the few mjy 
flux-density limit of our survey, enables the detection of a 
population of almost entirely starburst-dominated galaxies. 
This is also consistent with recent work by Trichas et al. 
(2009); although the 24/mi/X-ray initial selection in Trichas 
et al. (2009) implies that most objects in their sample host 
an energetically important AGN, the additional 70/im de- 
tection criterion results in these systems also being strongly 
star- forming. 
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